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We study the time evolution and gravitational wave emission properties of a black hole orbiting 
inside an accretion disk surrounding a massive black hole. We simultaneously solve the structure 
equations of the accretion disk in the presence of heating, cooling, and viscosity as well as the 
equations governing the companion orbit. The deviation from the Keplerian distribution of the 
angular momentum of the disk due to pressure and advection effects causes a significant exchange of 
angular momentum between the disk and the companion. This significantly affects the gravitational 
wave emission properties from the binary system. We show that when the companion is light, the 
effect is extremely important and must be taken into account while interpreting gravitational wave 
signals from such systems. 
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I. INTRODUCTION 

The study of gravitational wave emission from binary systems has received a significant boost in recent years because 
of the realization that the detection of gravitational waves would directly identify compact and strongly gravitating 
bodies, such as neutron stars and black holes. The Laser Interferometric Gravitational Wave Observatory (LIGO) and 
Laser Interferometer Space Antenna (LISA) project instruments are being constructed to achieve these goals In 
order to be able to obtain as accurate information about the radiating compact bodies as possible, efforts are being 
made to obtain correct forms of quadrupole radiation from a binary system In binary systems composed of 

only neutron stars and stellar black holes, these computations are adequate. However, when studying effects around 
a massive black hole which is assumed to be present in centers of many galaxies, one needs to consider an additional 
effect - the effect of an accretion disk. It is widely believed that galactic centers are endowed with massive black holes 
, and in order to explain the observed luminosity from a galactic core, one needs to supply matters ranging from 
a few hundredths to a few solar masses per year, some of which may be in the form of stars ||. Some of the stars 
could be compact, namely, neutron stars and stellar mass black holes which orbit the massive ones at the same time 
gradually spiraling in towards the center due to loss of angular momentum by gravitational waves. 

Chakrabarti |l(]] pointed out that the accretion disks close to the black hole need not be Keplerian and it would 
affect the gravitational wave properties. The radiation pressure dominated disks are likely to be super-Keplerian which 
would transfer angular momentum to the orbiting companion and in some extreme situations, can even stabilize its 
orbit from coalescing any further. This was later verified by time-dependent numerical simulations Jllj . When 
one considers the more general solutions of viscous, transonic, accretion disks ||l^,^3|, one finds that the angular 
momentum distribution close to the black hole could be sub-Keplerian as well, depending upon the viscosity and 
the angular momentum at the inner boundary of the disk. The disk becomes Keplerian roughly in a distance of 
XKep ~ (^ir) 2 from the black hole, where, M — v/a is the Mach number of the flow, v and a are the radial and 
sound velocities, and a 1 is a constant describing the viscosity [jl4|. Assume that a companion of mass is in an 
instantaneous circular Keplerian orbit around a central black hole of mass Mi. This assumption is justified, especially 
when the orbital radius is larger than a few Schwarzchild radius where the energy loss per orbit is very negligible 
compared to the binding energy of the orbit. The rate of loss of energy dE/dt in this binary system with an orbital 
period P (in hours) is given by |T^ , [l^1 

f = 3!< 1033 <iife ,s, W )1,a( ^ r '° /3eigs ""' (1) 
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where, 



MiM 2 
[l ~ Mi + M 2 

and 

M toi = Mi+M 2 . 
The orbital angular momentum loss rate would be, 



_dL 1 dE 

R aw - ~^\ gw - ^ (2) 



where 57 = ^/ GM\ / r 3 is the Keplerian angular velocity of the secondary black hole with mean orbiting radius r. The 
subscript 'gw' signifies that the rate is due to gravitational wave emission. In presence of an accretion disk co-planer 
with the orbiting companion, matter from the disk [with local specific angular momentum l{r)\ will be accreted onto 
the companion at a rate close to its Bondi accretion rate p7|Jlq] , 



M2 ~~ 7~2 i 2YU2 ( 6 > 

K V rel + a ) ' 

where p is the density of disk matter, A is a constant of order unity (which we choose it to be 1/2 for the rest of the 
paper), and v re i = Vdisk — VKep is the relative velocity of matter between the disk and the orbiting companion. The 
rate at which the angular momentum of the companion will be changed due to Bondi accretion will be |1C[ | , 

Rdisk = ~irUisk = M 2 (l Kep (x) - Idisk(x)) (4) 

Here, Ixep and Idisk are the local Keplerian and disk angular momenta respectively. The subscript on the left hand 
side signifies that the effect is due to the disk. If some region of the disk is sub-Keplerian {Idisk < liiep), the effect of 
the disk would be to reduce the angular momentum of the companion further and hasten coalescence. If some region 
of the disk is super-Keplerian, the companion will gain angular momentum due to accretion, and the coalescence is 
slowed down. In a thin disk with a high accretion rate, the Bondi accretion rate could be very high and the latter 
effect could, in principle, stop the coalescence completely JlO| |. 

In order to appreciate the effect due to intervention of the disk, we consider a special case where, M 2 << Mi and 
idisk << Ixep- ln this case, p ~ M 2 and M tot ~ M\. The ratio R of these two rates is, 

R = 5g± = 1.518 x l - 7 -^_a: 4 M 8 2 (5) 

Rgw Tio 

Here, x is the companion orbit radius in units of the Schwarzschild radius of the primary, Ms is in units of 10 8 Mq, 
pio is the density in units of 10 -10 gm cm -3 and Tio is the temperature of the disk in units of 10 10 K. It is clear that, 
for instance, at x = 10, and Mg = 10, the ratio R ~ 0.015 suggesting that the effect of the disk could be a significant 
correction term to the general relativistic loss of angular momentum. In the above example, both the disk and the 
gravitational wave work in the same direction in reducing the angular momentum of the secondary. Alternatively, 
when Idisk > 1-Kep they act in opposite direction and may slow down the loss of angular momentum [|To|] . In either 
case, the ratio R is independent of the mass of the companion black hole, as long as M 2 << M\. 

In what follows, we present equations governing the disk and the companion (§11). In §111, we solve these equations 
simultaneously along with eqs. (2) and (4) in a few typical cases which show varied nature of the disk structure and 
evolution of the companion. These disks are the generalization of the viscous, isothermal disks obtained earlier p2[ . 
We also present some interesting observations on non-axisymmetric disks containing spiral shocks. Finally, in §IV, we 
make concluding remarks. 



II. GOVERNING EQUATIONS 

To simplify the equations governing the accretion disks, we make the usual assumption that the disk is thin, so 
that the vertical averaging of density, pressure, and viscous stress, could be done and the vertical velocity component 
can be ignored. We assume the disk to be axisymmetric, an assumption relaxed in studying non-axisymmetric disks. 
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Instead of solving fully general relativistic equations, we assume the Paczyhski-Wiita potential approach which is 
sufficiently accurate to describe physical quantities around a Schwarzschild black hole, provided one is not too close to 
the horizon (r £ r g = 2GM\ /c 2 is the Schwarzschild radius of the primary black hole) . In this approach, the potential 
of the central body <&(r) = —GM\jr is replaced by a pseudo-Newtonian potential <I>pAr(r) = —GM\j(r — r g ). The 
effect we are presently discussing, namely the transport of angular momentum from the disk to the companion, is not 
of general-relativistic origin. Hence the major conclusions should not be affected by our simplified approach. 

The steady state accretion disk equations are: 
(a) The radial momentum equation: 



dx p dp x 3 



(6a) 



(b) The continuity equation: 



dx 



(pxhv) = 0, 



(66) 



(c) The azimuthal momentum equation: 



^ dljx) I d aPx 3 h dQ _ Q 
dx pxh dx &Kep dx 



(6c) 



(d) The entropy equation: 



ds 

EvT— = Q+ -Q-. (6d) 
dx 

Here Ixep and ^Kep are the Keplerian angular momentum and Kcplcrian angular velocity, respectively, £ is the 
density p vertically integrated, h — h(x) is the half-thickness of the disk at radial distance x, v is the radial velocity, s 
is the entropy density of the flow and Q + and Q~ are the heat gained and lost by the flow. We compute h(x) assuming 
the disk is in a hydrostatic balance equation in thethe vertical direction. l(x) is the angular momentum distribution 
of the disk matter. Here, we have chosen geometric units, thus x — r/r g is distance in units of a Schwarzschild radius, 
l(x) is in units of 2GMi/c, and velocities are in units of the velocity of light. We have implicitly assumed M 2 « Mi 
so that the gravitational effects due to the companion in shaping the disk could be ignored. However, locally, the 
companion is capable of exerting its effect to accrete matter from the disk, a £ 1 in the above equation is the viscosity 
parameter of Shakura and Sunyaev Jf4|, which is widely used to describe the viscous stress: w r< p = —aP. This stress 
transports angular momentum from the inner to the outer region of the disk. We choose total pressure (thermal plus 
ram) in this prescription in order that the angular momentum remains continuous across shock waves as well [L3| . 

The equation governing the companion, treated as a test particle in the field of the massive black hole is simply 

{ dt> x-1 x 2 [ ' 

In the following section, we present simultaneous solutions of the four sets of equations: (2), (4), 6(a-d) and (7) 
obtained by a very accurate fourth-order Runge-Kutta method [M . 

We note from the estimate of the ratio R in (5) that it is independent of the mass of the companion. However, the 
evolution time scale L / (dL / dt) of the companion orbit depends inversely upon its mass M 2 . Thus, even when the 
effect is very small R « I, the smaller companion will evolve so slowly that the number of cycles will be significantly 
affected. While integrating the above equations, however, we had to consider the computing ability of our machines. 
This constrained us to study cases with faster evolution time scale: we chose M\ = 1O 8 M0 and M 2 = IO 6 M0. Second, 
we note from (5) that the effect is directly proportional to the density of the gas in the disk, which in turn depends 
upon the accretion rate Mi of the primary. It is customary to express accretion rates in astrophysics in units of 
Eddington rates MEdd — ^GMim v j <jt ~ O.2Af Afg yr _1 , where m p is the proton mass and <jt is the Thomson 
cross section. Unlike the accretion process onto ordinary stars radiating from its surface, a black hole accretion 
process need not be limited by its Eddington rate. Since at the most r\ — 0.06 fraction of the rest mass energy is 
released by accreting matter on a Schwarzschild black hole a critical rate of M crit ~ MEdd/v ~ 16-Msdd is ver Y 
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reasonable to choose. For concreteness, we assume Mi ~ 70M cr it (i.e., M\ ~ lOOOAfgdd)- The effects we describe will 
be proportionately weaker when smaller rates are used. 

To keep the problem simple enough, we have considered the companion orbits to be circular. The general elliptic 
orbit is easily studied by including the evolution of the azimuthal coordinate in conjunction with eq. (7). This will be 
done later. Second, we assume that the companion orbit is co-planer with the disk, so that the companion is always 
immersed inside the disk. When it is not so, one has to include the fraction of time the companion is exchanging 
angular momentum with the disk and the effect would be proportionately reduced as well. In the case of lighter 
companions, the time scale of evolution is very long, and it is not unlikely to imagine that the orbits originally away 
from the equatorial plane will gradually lose the momentum component normal to the disk by repeated interaction 
pl| and eventually come to the plane of the disk much before our effects become important. In the case of massive 
companions, the evolution due to gravitational waves could be very rapid, and they may remain off the plane, as is 
possibly the case for OJ287 |||. 



III. SIMULTANEOUS SOLUTIONS OF THE GOVERNING EQUATIONS 

Before we present the results of our investigation, we discuss briefly what type of accretion disk solutions are 
expected. From (6c) (also see Chakrabarti ]l2|]) one observes that a weakly viscous disk starts deviating from Keplerian 
distribution very far away, whereas the strongly viscous disk remains Keplerian until close to the black hole. The disk 
cannot remain Keplerian very close to the black hole, as the velocity increases and the advection term [first term in 
(6a)] becomes important. Similarly, when the accretion rate is high or very low, the radiation pressure [ |19| or the gas 
pressure [ p3[ becomes important and the pressure term [second term in (6b)] cannot be ignored. Both of these terms 
were ignored in the study of Keplerian disks [|l^,^4). Also, in a Keplerian disk, angular momentum l{r) — Ixep is 
used independent of viscosity prescription. But generally, for some ranges of viscosity and accretion rates, this need 
not be true. Since our effects are non-zero only in non-Keplerian disks (eq. 4), it is essential that we include these 
effects. Near x = XKep where l{x) ~ Ixep, the distribution rapidly deviates from Keplerian to highly sub-Keplerian. 
This causes a 'micro-burst' of the gravitational wave emission as we shall show below. 

An important class of stable solutions of 6(a-d) involve shock waves |l2| , ^3| where the centrifugal barrier of the flow 
brakes the radial motion of the disk before the disk can continue through the sonic point to become supersonic, thus 
satisfying the boundary condition on the horizon. At the shock wave, the density, velocity and temperature change 
discontinously and the effect we are considering is expected to be discontinuous as well. This causes a 'micro-glitch' 
in the gravitational wave. In the case of non-axisymmetric disks, the spiral shocks cause micro-glitches to appear 
repeatedly depending on the number of spiral arms. This will be discussed towards the end of this Section. 

A few cases of the solutions are presented here, which cover all possible types of solutions. All the disk solutions are 
characterized by three parameters (instead of four, since the flow passes through one or more sonic points jl2| ) . These 
parameters are U n [angular momentum at the inner edge of the disk; the integration constant of equation (6c)], Mi 
(accretion rate on the primary; the integration constant of eq. 6b) and Xi n (the location of the inner sonic point; this 
defines specific energy of the flow at a given point). Alternatively, we could choose XKepi t ne location where the disk 
starts deviating from the Keplerian distribution, but we prefer to choose the sonic point location for convenience. We 
also choose an a, the unknown viscosity parameter, the polytropic index of the gas 7 which defines the specific internal 
energy of an ideal gas: e = (7— l)~ 1 P/p and the Q + — Q_, the relative importance of cooling and heating [eq. (6d)]. 
These are not completely independent parameters, but to obtain them one requires to include other equations in the 
list 6(a-d) to describe the viscous mechanism (such as poorly understood turbulence and convections), and cooling 
processes (such as the Compton effect, bremsstrahlung, pair creations and annihilations etc.). Instead of bringing in 
these equations we choose reasonable values for these quantities. Since the ratio R [eq. (5)] does not depend on M 2 
(but the orbital evolution time does) we consider only the case of Mi = 1O 8 M and M2 = 1O 6 M0 in order to hasten 
the evolution of the companion orbit. 

Case A: Figures l(a)-c) show results where the disk always remains sub-Keplerian after deviating from the Keplerian 
disk at XKep ~ 90r g . The disk smoothly passes through the inner sonic point at xi n = 2.3. Other parameters are 
kn = 1.7, Mi = lOOOMEdd, 7 = 5/3, a = 0.02, Q_ = Q + . In Figure la, we notice that the flow quickly becomes 
highly sub-Keplerian first. However, before entering the black hole it becomes only moderately sub-Keplerian. Figure 
lb shows the ratio R = Rdisk/ R gw - The ratio R jumps to almost 0.1 around x = 80 before decreasing to a very 
small value close to a black hole. Figure lc shows the number of times the companion orbits the primary (twice the 
number of full gravitational waves emitted). The solid curve is drawn including the effect of the accretion disk, while 
the dashed curve is drawn considering the usual binary orbit evolution (eq. 2) without the presence of the disk. Time 
in the abscissa denotes time passed since the companion entered the sub-Keplerian disk. Two effects are clear: (a) 
the binary coalescence takes place roughly 10% times faster and (b) the number of orbital cycles is also about 10% 
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times higher at the time of coalescence. If an accretion rate of M\ = MEdd were chosen instead, the effect would be 
reduced by a factor of 1000. If a lighter black hole of M2 ~ M Q was chosen instead, a longer orbital evolution time 
due to weaker gravitational wave loss gives rise to the same effect. 

Case B: In this case we choose a solution with a standing shock wave. The disk parameters are: U n — 1.6, Xi n — 2.87, 
Mi = lOOOAlEdd, 7 = 4/3, a = 0.05 and Q- = 0.5Q+- Figure 2a shows the Mach number variation as a function 
of distance from the black hole. The arrowed curve is followed by the disk after it deviates from Keplerian disk at 
XKep = 480r g . The disk first passes through the outer sonic point (located at x ou t — 50), then through the shock 
at x s = 13.9 and finally enters the black hole through the inner sonic point at Xj„ = 2.87. The shock location or 
the location of the outer sonic point is not a free-parameter, but are self-consistently determined from the Rankine- 
Hugoniot relation (2^Jl^]. The shock solution is always chosen if it is available to the disk, since the entropy at the 
inner sonic point is higher compared to its value at the outer sonic point, and the required entropy must be generated 
at the shock. Figure 2b shows the ratio of disk to Keplerian angular momentum distributions. Figure 2c shows the 
ratio R as a function of the distance. The ratio becomes almost 5, (5 x 10~ 3 for Mi = MEdd) a micro-burst of a sort, 
around x = 400r s . There is also a glitch at the shock location. In cases with a stronger shock wave the glitch would 
be stronger. 

Case C: In this case we choose disk parameters so as to obtain a super-Keplerian region in the disk. We choose 
l in = 1.88, x c = 2.2, Mi = l000M E dd, a = 0.005, 7 = 4/3 and Q_ = Q+. The disk deviated from Keplerian disk at 
XKep = 7.5r g . Figure 3a shows the ratio of disk to Keplerian distributions which clearly shows the sub-Keplerian as 
well as super-Kcplerian regions. Figure 3b shows the ratio R (eq. 5). The fractional change in orbital cycle number 
with and without the disk is ^ ~ R. Thus, SN ~ 1 only when N ~ R^ 1 ~ | dt/dt^ • Thus, even if R is small, lighter 
companions should survive long enough to feel the effect of angular momentum exchange. Chakrabarti |l0[ ] considered 
a thin disk where the density of the disk was high enough to stabilize the companion orbit in the super-Keplerian 
region. 

Case D: In this case we solve non-axisymmetric disk equations |^6| instead of eqs. 6(a-d). Here spiral shocks formed 
would produce repeated glitches in the gravitational waves pattern. The simplest solutions of the non-axisymmetric 
shocks are obtained by assuming self-similarity in x and all the disk velocity components vary as q^cfyx^ 1 ^ 2 and the 
density of the disk varies as q p ((j))x~ 3 ^ 2 . Here, azimuthally varying coefficients q%{<f) and q p {<f>) are to be determined 
from boundary conditions. Figure 4a shows a typical solution for the velocity coefficients when a two armed spiral 
shock solution is considered. The solid, long-dashed and short-dashed curves show the the radial, azimuthal and 
sound velocity coefficients and the dotted curves show the density coefficient as they vary with the azimuthal angle. 
The shocks are located at <p = and <j) = n. In a Keplerian disk, the azimuthal velocity coefficient would be unity 
throughout the disk. In this example, the azimuthal velocity coefficient varies from 92% of Keplerian to 26% Keplerian 
as the flow crosses the shock front. Other components also suffer a jump. In a single circular orbit, the companion thus 
passes twice through these jumps. Figure 4b shows (in arbitrary units) the glitches in the ratio R in a single orbit. In 
an axisymmetric disk, the glitch appears only once, but in a disk with spiral shocks the effect occurs repeatedly and 
cumulative effect becomes important due to the repeated passage of the companion through the shock. This could 
significantly modify the shape of the gravity wave signals. 



IV. CONCLUDING REMARKS 



It is widely recognized that accurate templates of possible signals may be essential to determine the nature of 
radiating compact bodies Q]. In this paper, we have discussed several important ways a gravitational wave signal 
from a binary companion could be modified in presence of an accretion disk. We find that even under very normal cir- 
cumstances, the effects will be sufficiently significant and our effect may influence the templates constructed assuming 
the absence of accretion disks. 

In a binary system containing lighter mass black hole components, the accretion disk need not be present. Systems 
involving massive black holes at the galactic center may necessarily contain accretion disks and lighter black hole 
companions. The frequency of the gravitational wave f gw = 2.25 x 10" 4 x~ 3/2 M 8 _1 is well outside LIGO sensitivity, 
but could be well within LISA sensitivity j|,^7]]. By self-consistently solving the equations governing the accretion 
disk structure and the evolution of the binary orbit, we first showed that accretion disks close to the black hole are 
in general non- Keplerian. In the sub-Keplerian region of the disk, the residence time of the companion inside a disk 
and the probability of its observation would be reduced. On the other hand, the super-Keplerian region enhances the 
residence time and the probability of observing these systems is higher. We also find that the orbital evolution may 
be faster away from the black holes where the disk angular momentum distribution starts deviating from Keplerian 
distribution. The population density of compact stars close to galactic nuclei should be affected by their interaction 
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with the disk. These effects shouid be taken into consideration while determining the band of maximum sensitivity 
of future instruments for gravitational wave astronomy. 

The discussions made in this paper involving a companion black hole are valid even when a neutron star is chosen 
instead and may approximately remain valid when an ordinary star orbits the central black hole. In the latter case, 
the angular velocity of the disk changes significantly along a radial direction across the star. This would cause 
some angular momentum of the disk to spin up or spin down the star itself rather than changing its orbital angular 
momentum. Furthermore, the star may lose angular momentum through winds. Therefore, our result need not be 
strictly valid in these systems. These effects are negligible if the companion is a black hole or a neutron star because 
of its small size and the absence of winds. 

The assumption of a thin disk in vertical equilibrium (namely, that the vertical velocity is negligle compared to the 
radial or azimuthal velocity) enabled us to integrate the governing equations. Numerical simulations of fully three 
dimensional disks p8[ indicate that the assumption of the vertical equilibrium adequately describe the disk properties 
and therefore we do not believe that the conclusions drawn in the present paper are affected by this assumption. 
Another implicit assumption has been that the disk remains continuous (and does not break apart in the form of 
rings as in the case of orbiting matter around Saturn) even in the presence of an orbiting companion. The formation 
of gaps in the disk is possible only if the instantaneous gap is not filled in by the accreting matter through radial 
pressure or viscous forces 21,23]. This implies that either the Roche radius (Rl) of the star orbiting at radius r is 



greater than the disk thickness, Rl ~ (Afe/Mi) 1 /^ £ h ~ ar 3 / 2 ~ 0.5r (here, the sound speed a ~ l/y/Sr^ 1 / 2 ), or, 
the viscosity parameter is so small that the angular momentum transfer rate by tidal coupling through the satellite 
is higher than that by viscosity: a < l/40(M 2 /M 1 ) 2 (r/h) 5 ~ 3/8(M 2 /M 1 ) 2 . It is clear that for the cases we are 
interested in, namely, for M2/M1 ^ 1CU 6 , and a £ 1CU 3 , neither of these conditions would be satisfied. Thus, we do 
not think that gaps would be formed by orbiting black holes or neutron stars. 

Our goal in this paper has been to indicate a new physical effect which may change the gravitational wave pattern 
significantly. Construction of accurate templates for inferring component masses of the gravitating systems is beyond 
the scope of the present paper. It is possible that one could estimate the mass of the central black hole by comparing 
the observed optical and UV radiation spectra and the hard and soft X-rays with the theoretically derived spectra using 
these general disk models 13 3(]]. The accretion rate of the Keplerian disk could be inferred from the normalization 



of the optical and UV flux as well. Thus two parameters are easily eliminated. Viscosity is not an well understood 
process in the context of accretion phenomena, but typical values of the parameters have been presented in the 
literature from time to time |^l| , ^2| which we have considered here for simplicity. Since whether the flow is Keplerian 
or non-Keplerian depends very crucially upon the viscosity parameter ]l3[ | , undoubtedly, a complete resolution of the 
present problem is hinges upon a better understanding of the viscosity of the accretion disk. 
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FIGURE CAPTIONS 



FIG. la: Ratio of disk angular momentum distribution to Keplerian distribution of a disk which is entirely sub- 
Keplerian for x < XKep — 90. See text for flow parameters. 

FIG. lb: Ratio of the rates of change of angular momentum of the companion due to exchange with the disk and due 
to gravitational wave emission. The ratio is highest in regions closer to the Keplerian boundary. 

FIG. lc: Comparison of number of orbital cycles in a binary with disk (solid) and without disk (dashed) as time passes 
since the companion enters the sub-Keplerian region of the disk. Companion falls faster when the disk is present. 



FIG. 2a: Variation of Mach number of the disk which includes a shock wave at x s ~ 13.9. The arrowed curves are 
the solutions chosen by the flow. See text for flow parameters. 

FIG. 2b: Ratio of disk angular momentum distribution to Keplerian distribution of the disk with a shock which is 
entirely sub-Keplerian for x < XKep = 480. 

FIG. 2c: Ratio of the rates of change of angular momentum of the companion due to exchange with the disk and due 
to gravitational wave emission. The ratio is highest in regions closer to the Keplerian boundary. Note the glitch at 
the shock location which could be very high for stronger shocks. 

FIG. 3a: Ratio of disk angular momentum distribution to Keplerian distribution of the disk where the sub-Keplerian 
disk below x < XKep — 7.5 becomes super-Keplerian close to the black hole. See text for flow parameters. 

FIG. 3b: Ratio of the rates of change of angular momentum of the companion due to exchange with the disk and due 
to gravitational wave emission. The ratio is highest in regions closer to the Keplerian boundary. Note the change in 
sign of the ratio as the companion enters the super-Keplerian region. 

FIG. 4a: Velocity and density variations with azimuthal angle in a non-axisymmetric disk with a two-armed spiral 
shock waves located at <f) — and <f> = n. The velocities are in units of local Keplerian velocity, while the density is 
in an arbitrary unit. 

FIG. 4b:Ratio of the rates of change of angular momentum of the companion due to exchange with the disk and due 
to gravitational wave emission. The jump in the ratio at the spiral shocks produces glitches twice per orbital cycle 
(once per gravitational wave signal). 
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